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We study collective electronic excitations in graphene in the integer quantum Hall regime, concen- 
trating mainly on excitations with spin reversal such as spin-flip and spin-wave excitations. We show 
that these excitations are correctly accounted for in the time-dependent Hartree-Fock and strong 
magnetic field approximations, in contrast to spin-conserving (magneto-exciton) modes which in- 
volve a strong Landau-level mixing at non-zero wave vectors. The collective excitations are discussed 
in view of prominent theorems, such as Kohn's and Larmor's. Whereas the latter remains valid in 
graphene and yields insight into the understanding of spin-dependent modes, Kohn's theorem does 
not apply to relativistic electrons in graphene. We finally calculate the exchange correction to the 
chemical potential in the weak magnetic field limit. 

PACS numbers: 78.30.Na, 73.43.Lp, 81.05. Uw 



I. INTRODUCTION 

The role of electron-electron interactions in graphene 
(two-dimensional graphite) is still a debated issue. 
Whereas most of its electronic properties can be un- 
derstood within a model of two-dimensional (2D) non- 
interacting massless Dirac fermionsP there are some 
experimental indications for the presence of Coulomb 
interactions.^ These correlations, which may be quan- 
tified by the graphene fine-structure constant aa = 
e 2 /hev-p ~ 2.2/e, in terms of the Fermi velocity vp and 
the dielectric constant e, seem, however, to be weak and 
long-ranged. Therefore, strongly-correlated phases that 
are expected in the large-ac limilP or for short-range 
Hubbard interactions^^ are unlikely to occur in undoped 
or moderately doped graphene. Theoretically, a pertur- 
bative Fermi-liquid-type treatment of the Coulomb in- 
teractions yields a logarithmic divergence of the Fermi 
velocity^ a renormalization of thermodynamic quantities 
such as the compressibility^ as well as to a control of the 
orbital magnetic susceptibilityP 

The situation is different if the graphene electrons are 
exposed to a strong magnetic field that quantizes their 
kinetic energy into non-equidistant Landau levels (LLs), 
e„ = (Afi.WF/^s)v / 2n, where A = ±1 is the band index, 
Ib — \Jhcj eB is the magnetic length, and n denotes the 
LL indexpJ The most prominent consequence of this rela- 
tivistic LL quantization and the presence of a zero-energy 
LL for n = is a peculiar integer quantum Hall effect 
(QHE), with an unusual sequence of Hall plateaus! 10 * 11 1 
The LLs are highly degenerate, as in the usual 2D elec- 
tron gas (2DEG) with a parabolic band dispersion, where 
the density of states per LL (and per unit area) is given 
by the flux density rig = l/2wl B = eB/h, which is pro- 
portional to the perpendicular magnetic field. The filling 
of the LLs is then characterized by the ratio (filling fac- 
tor) v = n e i/riB between the 2D electronic density n e i 
and Us- A partially filled LL may then be viewed as 



a strongly-correlated electron system with a quenched 
kinetic energy, and its most prominent manifestation is 
a fractional QHE that has recen tly been observed in 
suspended graphene samples! 12 * 13 ! Prior indications for 
strong interactions stemmed from the high-field QHE at 
v = 0, which indicates a stronger lifting of the four- fold 
spin- valley degeneracy of graphen e tha n what one would 
expect from single-particle effectsP^ESl 

Similarly to the B = case, the Coulomb inter- 
action between electrons in completely filled LLs may 
be viewed as a weak perturbation because of the en- 
ergy gap between adjacent LLs. Rs role in the dis- 
persion relation and Fermi velocity renormalization of 
graphene is an open question which has been addressed 
both theoretically-^"^ and experimentally, in the frame- 
work of transmission spectroscopy!^^!] As compared 
to the 2DEG with a parabolic band dispersion (as in 
GaAs/AlGaAs and Si/SiGe heterostructures), the situa- 
tion is strikingly different in graphene, where the effect of 
electron-electron interactions may be probed at zero wave 
vector. Indeed, in the former, LL quantization leads to 
a set of equidistant LLs separated by the cyclotron fre- 
quency. Kohn's theorem states that in these systems, 
homogeneous electromagnetic radiation can only couple 
to the center-of-mass coordinate. Therefore internal de- 
grees of freedom associated with the Coulomb interaction 
cannot be excited by such optical probesP^Then, the dis- 
persion relation of spin-conserving magnetoplasmons at 
zero wave-vector is equal to the bare cyclotron energy, 
irrespective of existing electronic correlations.^ A simi- 
lar consideration holds also for spin wave (SW) modes, 
for which Larmor's theorem states that the Coulomb in- 
teraction does not renormalize the zero-wave-vector dis- 
persion of the spin excitonsP^ However, the dispersion 
of spin-flip (SF) modes in a 2DEG are shifted from the 
bare cyclotron resonance even at zero wave-vector, due 
to electron-electron interactions. Therefore, these excita- 
tions are the only suitable modes to study the many-body 
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effects in a 2DEG by means of optical measurements 

Furthermore, electron- electron interaction in the 
regime of the integer QHE yield collective excitations 
that are different from those in the 2DEG - instead of 
inter-LL excitations with a rather weak wave-vector dis- 
persion, called magneto-excitons (ME), one finds linear 
magneto-plasmons that involve superpositions of differ- 
ent LL transitions! 30 ^ Also MEs that may play a role 
in the vicinity of q = have been studied theoretically 
in graphene, and it has been shown that Coulomb in- 
teractions yield a reno rmalizat ion of the transition en- 
ergy at zero wave vectoJ^ 17 * 32 ! that indicate that Kohn's 
theorem does not apply to graphene. In comparison to 
these works, here we put more emphasis on spin-changing 
modes and on the effect of LL mixing. 

A convenient way of assessing the magnetic-field 
strength is in terms of four characteristic length scales: 
the magnetic length Ib, the carbon-carbon distance o, 
the Fermi wavelength Xf and the Thomas- Fermi screen- 
ing length Xtf ~ ^fZ&g where qq = e 2 /evp measures 
the relative strength of Coulomb interactions. In prac- 
tise, the magnetic length is always much larger than the 
lattice spacing lg ^> a because the flux per unit cell is 
much smaller than the flux quantum (B <C 40000 T). 
The weak LL mixing approximation - which we will use 
when studying the particle-hole excitations - requires 
e / (eIbojc) 1, where cue is the cyclotron frequency, 
and corresponds to a strong field such that lg <C Xtf- As 
in graphene ag is of order 1, Xtf ~ Xf and the weak LL 
mixing is also the small filling factor limit Ib <C Xf- In 
the following, we will consider two limits: either a strong 
magnetic field (meaning Ib <C Xp or typically B 3> 20 T) 
or a weak magnetic field (meaning Ib 3> Xp or typically 
B « 20 T). 

In this paper, we study both spin-conserving ME and 
spin-dependent SF and SW modes in the regime of the 
integer QHE. Following the scheme introduced by Kallin 
and Halperin (KH) for the 2DEGP 3 the Coulomb in- 
teraction is treated within the framework of the time- 
dependent Hartree-Fock (TDHF) and strong-field ap- 
proximation [for which e 2 /(eIb^c) 1; which insures 
that LL mixing is weak], the validity of which is dis- 
cussed. Indeed, we find that whereas SF and SW may be 
accounted for correctly in graphene within the KH ap- 
proximation in the limit of a strong magnetic field and 
when the Fermi level lies near the n = LL, its validity 
in the treatment of spin-conserving ME is questionable 
even in these limits. This difference between MEs and 
SF (and SW) excitations stems from the depolarization 
term, accounted for in the random-phase-approximation 
(RPA), which is present only in the ME dispersion and 
which yields a strong LL mixing at non-zero values of 
the wave vector. This LL mixing eventually leads to the 
formation of linearly dispersing plasmon-type modes that 
have been obtained within an R PA tre atment of the elec- 
tronic polarizability in graphene! 30 * 31 ! 

Finally, we consider the opposite limit of graphene in 
a weak magnetic field, and compute the exchange correc- 
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FIG. 1. Sketch of the particle-hole excitations studied in the 
text. Each Landau level is split in two sub-levels, separated by 
the Zeeman gap, AE Z = q^lbB. Here, f indicates s z — +1/2 
and \., s z — —1/2. Label ME stands for magneto-exciton (or 
magneto-plasmon) where the particle and the hole have the 
same spin. SF denotes the spin-flip excitation, in which the 
electron and the hole do not only reside in different LLs but 
also have different spin. Finally, SW denotes the spin-wave 
mode, which is an intra-LL transition where the electron and 
the hole have a different spin orientation and where we have 
taken into account a non-zero Zeeman energy. 



tion to the chemical potential. We find that the exchange 
correction to the single-particle dispersion presents the 
same dependence in the two limits, strong and weak mag- 
netic field, proportional to the square root of the ultra- 
violet cutoff in the Landau levels. On the other hand, 
the exchange correction to the particle-hole dispersion 
diverges logarithmically with the cutoff. 

The paper is organized as follows. In Sec. |TT1 we first 
revisit the Kohn's (Sec 



II A) and Larmor's (Sec. II B) 



theorems in the context of graphene and their impact 
on collective excitations in general. We then study the 
excitonic modes in graphene in a strong magnetic field, 
within the KH approximation. In Sec. |III| we calculate 
the exchange correction to the chemical potential, in the 
weak magnetic field limit and/or for highly doped sam- 
ples. Our main conclusions are summarized in Sec. |IV[ 
and the technical details of the calculations are provided 
in the appendices. 



II. EXCITONIC MODES IN GRAPHENE IN 
THE INTEGER QHE 

In spin-flip modes, an electron is both promoted from 
one LL to the next one and its spin is reversed. They 
carry a spin S z = s\ — s h z — ±1, where st is the z- 
component of the electron (hole) spin. We use the term 
magneto-exciton (ME, sometimes also called magneto- 
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plasmon in the literature) to denote spin-conserving ex- 
citations, where the electron and the hole reside in dif- 
ferent LLs and have the same spin. In SW modes, the 
two particles have the same LL index but opposite spin. 
SF modes, MEs and SWs are the basic excitations of a 
quantum Hall system with a finite Zeeman splitting, as 
sketched in Fig. [T] 

A. The fate of Kohn's theorem in graphene 

Before calculating the dispersion relations of the dif- 
ferent excitonic modes, we discuss here qualitatively the 
expectations for graphene with respect to the 2DEG 
with a par abolic band dispersion. In the latter, Kohn's 
theorem^ states that electromagnetic absorption, irre- 
spective of the strength of the Coulomb interaction, oc- 
curs only at the cyclotron frequency uc — eB/nif,, in 
terms of the band mass m^. Here and in the remain- 
der of this paper, we consider a system of units with 
h = c = 1. Because this absorption process is associ- 
ated with an inter-LL transition from the last occupied 
to the first unoccupied level, it means that the lowest- 
energy ME must converge to the non-interacting value 
at zero wave-vector. This statement remains valid also 
for higher harmonics, such that the ME dispersion has 
^m_b(<? — > 0) — > mu>Ci where m — n e — rih is the differ- 
ence between the LL index of the electron (n e ) and that 
of the hole (n/j). 



1. Kohn's theorem in the 2DEG 

In order to investigate the fate of Kohn's theorem 
in graphene, let us recall the main steps in the argu- 
ment for the 2DEG. In the latter, the Hamiltonian of 
N non-interacting electrons can be expressed in terms of 
the gauge-invariant momenta II j — pj + eA(rj), where 
rj and Pj are the position and its conjugate (gauge- 
dependent) momentum, respectively, of the j-th electron 
with charge — e (we choose e > to be the positive ele- 
mentary charge), 

3=1 

From the total gauge-invariant momentum II = IX, , 
we define the raising and lowering operators II ± = H x i 
illy, which satisfy the commutation relations [II , Ho] = 
=F(l/mbZ^)n ± with the Hamiltonian Ho. This leads to 
the equation 

H (ll ± \^)) = (E ±uc)(ll±\i> )), (2) 

which means that the application of on an (JV- 
particle) eigenstate \ip a ) (with energy E°) of H yields 
another eigenstate with energy E i ojc- 



The first observation is that this equation remains valid 
also in the presence of electron-electron interactions V 
that commute with the total momentum [II, V]=0, such 
as the Coulomb interaction, if one replaces the non- 
interacting state by an eigenstate of the full 
Hamiltonian H = Hq + V, as well as the energy E° by 
that, E, of the state \ip). 

Second, one notices that the electromagnetic light field 
with frequency lo couples to the electronic system via the 
Hamiltonian 

H LM (t) = ^e-^'EH + H ' c -' ( 3 ) 

3 

where E(u>) is the electric component of the light field and 
Vj the velocity operator of the j-th electron. In the 2DEG 
with a parabolic band dispersion, the velocity operator 
is readily expressed in terms of the gauge-invariant total 
momentum, 53j v i = H/mj, such that the light-matter 
coupling ([3]) is linear in the operators TV* 1 . As mentioned 
above, this induces then a transition from a state \tp) 
with energy £ to a state n~|^>) with energy E±uic, i.e. 
the only absorption peak for light occurs at the cyclotron 
frequency wc- 

2. Difference in graphene 

Although also in graphene the total gauge-invariant 
momentum II commutes with the interaction Hamilto- 
nian V but not with Hq, one first notices that it may 
no longer be expressed in terms of the velocity opera- 
tors of (now relativistic) electrons because of the vanish- 
ing band mass. The velocity operator is a 2 x 2 matrix 
Vj = vpa-j = VF(<7j,aj), in terms of the Pauli matrices 
a x and a v , and it is not a conserved quantity even in the 
absence of interactions. The application of the velocity 
operator on an eigenstate of the non-interacting Hamilto- 
nian H = J2 3 VF[Pj+eA(rj)]-<Tj, for which [vj,H ] ^ 0, 
yields, even in the absence of a magnetic field, sponta- 
neous inter-band transitions that are at the origin of the 
so-called zitterbewegungW^ As a consequence, the light- 
matter coupling Hamiltonian the form of which is 
also valid for graphene, may no longer be expressed in 
terms of n*. Indeed, the velocity operator in Hamilto- 
nian ([3| yields transitions involving LLs with adjacent 
indices n and nil, as in the 2DEG, but the zitterbewe- 
gung translated to the magnetic-field case^ furthermore 
yields inter-band excitations, such that the dipolar se- 
lection rules Xh , n — > A e , n i 1 are associated with the 
energies 

Etn°' Xh) = ^ (Ae^^ + T - \HVK) , (4) 

where one expects absorption peaks. Therefore, al- 
ready in the non- interacting limit, one expects a plethora 
of absorption peaks, that have indeed been observed 
experimentaHyp^H^l and not a single cyclotron resonance 
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as in the case of the 2DEG with a parabolic dispersion 
relation. 

Furthermore, because the kinetic Hamiltonian ([!]) be- 
comes Hq = v-p Y] j IIj • (Tj in graphene, one loses the 
possibility of writing an equation of the type ^ for 
graphene, neither in terms of the total momentum II nor 
with the help of Y]j Vj , which as we mentioned is not 
conserved. There is thus no protection of the energies 
Q when interactions are taken into account. Indeed, 
the latter renormalize the absorption energies] 16 * 17 * 32 ! as 
we discuss below, in contrast to the 2DEG, where the 
absorption energy is protected by Kohn's theorem, and 
the ME modes no longer converge to the non-interacting 
inter-LL transition energies 

E un nh) = ^V2(X e ^~ e -X h ^) (5) 

IB 



which all spins are spontaneously polarized. This fer- 
romagnetic state arises due to exchange-interaction ef- 
fects when not all subbranches of a particular LL are 
completely filled.^ The Goldstone mode is characterized 
by a dispersion relation that vanishes (as q 2 for a SW 
mode^) in the zero-wave-vector limit, ojg(q — > 0) — > 0, 
which means that the different states of the ground-state 
manifold (i.e. the different polarizations) are connected 
by a global rotation of zero energy cost that is precisely 
the q — Goldstone mode. In the presence of the Zee- 
man effect, which chooses a particular orientation of the 
total spin, one thus obtains a SW mode that tends to the 
energy f2sw(<7 —> 0) — » g^sBSz, where S z = 1, as stated 
by Larmor's theorem. 



in the zero- wave- vector limit. 



B. Larmor's theorem applied to graphene 

In addition to ME excitations that do not involve the 
electronic spin, one may investigate spin excitations on 
rather general grounds. Larmor's theorem states that in 
the long- wavelength limit, the SW dispersion tends to 
the (bare) Zeeman splitting, ilsw{q 0) — > /////;/>'. 2li 
This theorem may be understood from the symmetries 
of the Hamiltonian H = H + H int + Hz- In the absence 
of the Zeeman term Hz, the Hamiltonian respects the 
SU(2) symmetry associated with the electronic spin, i.e. 
both the total spin operator S 2 ot and any of the com- 
ponents S^ ot , for fi — x,y, z, commutes with the Hamil- 
tonian. Since one cannot diagonalize all components of 
the total spin simultaneously, one needs to choose a par- 
ticular one, and this is naturally the one chosen by the 
Zeeman effect (here S z ot ), such that the full Hamilto- 
nian commutes with Sf ot and Sf ot . The quantum num- 
bers associated with the spin, S and S z , are therefore 
good quantum numbers for the full interacting A^-particle 
Hamiltonian, such that all possible states have energies 
E = E(S, S z , ...) + g[iBBS z , where the dots ... represent 
other quantum numbers that characterize the interacting 
system. The essence of this expression is that the full in- 
teracting ./V-particle system may be viewed as a large spin 
that precesses in a magnetic field with the fundamental 
(Larmor) frequency ujl = g^BB. Whereas this frequency 
is affected by the (crystalline) environment via the effec- 
tive g-factor, the latter remains unaltered by the electron- 
electron interactions. Applied to the present problem of 
collective excitations, this means that the Zeeman term 
does not represent a further complication to the SU(2) 
symmetric Hamiltonian H + H int , which thus needs to 
be diagonalized first. 

These rather obvious considerations allow us to un- 
derstand easily Larmor's theorem if one notices that, in 
the absence of a Zeeman effect, the SW mode is just 
the Goldstone mode of a ferromagnetic ground state in 



One notices that, in contrast to the above discussion 
of Kohn's theorem, the (non-)relativistic character of Hq 
has never played a role in the argument, and Larmor's 
theorem therefore also applies in the case of graphene. 
Moreover, one is confronted in graphene with an ad- 
ditional two-fold valley degeneracy, that may be taken 
into account by an SU(2) valley isospin. Although the 
SU(2) valley symmetry is not respected by the inter- 
action Hamiltonian, the symmetry-breaking terms are 
strongly suppressed (by a factor of a/Is ~ 0.005-y/S[T], 
in terms of the carbon-carbon distance a — 0.14 nm) 
such that the interaction Hamiltonian is approximately 
SU(2) valley-symmetric! 37 * 38 ^ Therefore the above argu- 
ments apply also to possible valley-ferromagnetic states 
in graphene, i.e. there are valley-isospin-wave modes that 
vanish in the q — > limit and that may become eventu- 
ally gapped by a "valley-Zeeman" effect H v _z that, if 
it may be written in terms of components of the total 
valley-isospin, yields a simple energy offset to the disper- 
sion. In the remainder of the paper, we concentrate on 
collective excitations that involve only the physical spin. 



In addition to this generalization of Larmor's theorem 
to the valley isospin, it may also be generalized to the 
SF modes, which involve not only different spin states 
but also different LLs. The dispersion of the collective 
SF modes may be fully understood from the Hamilto- 
nian Hq + Hi n t, whereas the energy g^,BBS z associated 
with the Zeeman effect can be simply added at the end 
of the calculation as a global (wave-vector independent) 
constant. However, as we shall discuss below, the en- 
ergy of the SF modes does not converge to the simple 
sum of the Zeeman and the transition energies (|5| in the 
zero- wave- vector limit, but they are renormalized by the 
interaction energy, both in graphene and in the 2DEG. 27 
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C. Dispersion relation of the excitonic modes 

In graphene, the energies of ME, SW and SF modes 
can be expressed as: 



n ME ( q ) = Et' nh) +AE^^\ q ) 
^sw(q) = 9^bBS z + AE^^ n '^\q) 

9HB 



n SF (q) = +gy. B BS z + AE^^\q) 



(6) 
(7) 

) 

(8) 




A^=0;v5=l,v?=0 



f. \ [39] jj ere j?„ = Zj. ^2 - g xj^ FIG. 2. Dispersions (in units of e /sIb) of the excitonic modes 

= 1 and ^ J = 0. SW 



where 5 Z is the z-component of the exciton spin, and 
^kin ^ S ^ ne transition energy in the absence of interac- 
tions given by Eq. ((H). The contribution AE < - n °< s *'< nh < s i') 
consists of three terms (see Appendix[A|for details): a de- 
polarization or exchange term E x (q) , which is accounted 
for in the RPA approximation, a direct Coulomb in- 
teraction between the electron and hole (vertex correc- 
tions) E v (q), and the difference between the exchange 
self-energy of the electron and that of the hole, E exc h = 
E e — S/j. Notice that E x (q) is only relevant for the ME, 
because only particles with the same spin can be recom- 
bined by means of electron-electron interactions. 

It must be kept in mind that, in a 2DEG, the RPA 
term, which determines the maximum of the ME disper- 
sion at a wave- vector q ~ 1/-Rc m the TDHF approxi- 
mation, mixes different LLs, with a mixing amplitude on 
the order of e 

clotron radius, where the Fermi momentum in terms of 
the index N F of the topmost fully occupied LL is k F — 
\j2Np + l/l B for a 2DEG and k F = ^/2N F + 5 NFfi /l B 
for graphene. This needs to be distinguished from the 
LL mixing at q = 0, which determines the stability of 
the LLs in the presence of electron-electron interactions 
and which scales as e 2 / eRc^c ■ Although the stability 
of the LLs in graphene is determined by the ratio be- 
tween the Coulomb energy e 2 /eRc and the LL separa- 
tion A„ = (V2v F /l B )(VN F + 1 - v^Vy7) ~ vp/Rc^ 
which happens to be the scale-invariant fine-structure 
constant olq = 6 2 /svf, the situation is again different 
at the maximum of the ME dispersion at q ~ I/Rq- 
The order of magnitude of the q =/= LL mixing in 
graphene may be obtained by replacing in e 2 / (eIb^c) the 
2DEG cyclotron frequency uoc — eB/rrib by the density- 
dependent cyclotron frequency u>c{n) — gBv 2 f / /j,, where 
/i = (vf/Ib)\/2N f is the chemical potential. As a con- 
sequence, the validity of the KH approximation fails not 
only for weak magnetic fields, as in the standard 2DEG, 
but also at high and intermediate filling factors because 
the effective cyclotron frequency in graphene decreases 
as the number of filled LLs increases, leading to an in- 
crease of the LL mixing. Therefore, strictly speaking, the 
results of this section will be valid only in the strong- B 
limit and for N F near 0. However, we will see that the 
KH approximation can still be applied for spin-dependent 
excitations (SW and SF) slightly away from half-filling, 
but not to spin- conserving modes (ME) . This is a conse- 
quence of the absence of the depolarization term E x (q), 



4 

(a) 



n= 1 



n = 0™T"SW 



n = -1 




00 



studied for v — 0, 



N F = 0, 



(dotted green line), MEi,2 (dashed blue and red lines, re- 
spectively) and SFi,2 (solid blue and red lines, respectively) 
are represented. The thin horizontal line represents the dif- 
ference in kinetic energy between the electron and the hole 
Ekin ■ We have used for the Zeeman term an unphysically 
large value Q^bB — (l/10)(e 2 /eIb), for illustration reasons, 
(b) Schematic representation of the excitonic modes studied. 
Notice that MEi and ME2 are degenerate in the N c — > 00 
limit, as well as the SFi and SF2. 



which is the main source of LL mixing, in SW and SF 
modes, whereas it constitutes the main contribution to 
the dispersion of ME modes. 

After these general considerations on collective exci- 
tations, we now turn to a discussion of the modes at 
particular integer filling factors, which are described by 



W F - 2 + 2 (v. 



t 



N F 



(9) 



where N F is the index of the top most fully occupied LL, 
< < 1 is the filling of the spin-u branch of the n-th 
LL, and the factor of 2 accounts for the two-fold valley 
degeneracy of each spin branch. 



D. Modes at filling v = 



() 



At the charge neutrality point (for a filling factor v = 
0), the Fermi level is in the n = LL (i.e. Np — 0), with 
the spin-j. branch completely filled (z/P = 1) and an empty 
spin-f (z/£ = 0). The dispersion of the excitonic modes for 



this situation is shown in Fig. 2(a) The transitions cor- 
responding to the different excitations are schematized 
in Fig. 2(b)| To more easily distinguish between the 
different modes, we use the notation AE^ \ w n„(q). 

Therefore, the dispersion of the magnetoexciton modes 
ME 12 , Eq. ^ will correspond to the kinetic particle- 
hole energy difference plus a renormalization due to 
electron-electron interactions, AEq 1 ^' 2 (q), which reads 



J 0;l,0 ' 



A £(l,-l/2;0,-l/2)( g ) 

MEi and AE( Q - +1 / 2 ^ +1 ^ {q) = FJ^ 2 + V ( 



-^o,i.o(?)+4^o:i,o(?) f or 



0,-l;0,-l 



(?)H 



0,-l;0, 



_i(q) for ME 2 , where the expressions for Eq.j 



are given in Appendix|Bj Here V x (q) are matrix elements 
of the Hartree term, in which a particle-hole pair re- 
combines, exciting a new particle-hole pair (the usual 
bubble diagrams). On the other hand, V (q) is the 
Fock term, which accounts for the direct interaction of 
the excited electron and hole (ladder diagrams). No- 
tice that, due to particle-hole symmetry at this filling, 
AE^iq) = AE^(q) and the two modes are degen- 
erate. The first thing one notices is that the dispersion at 
q = is s hifted with respect to E^' [horizontal line in 



Fig. 2(a) . This is a consequence of the non-applicability 
of Kohn's theorem in graphene, as discussed in Sec. |II A| 
whereas in the 2DEG the theorem is satisfied due to 
a cancelation between the exchange self-energy and the 
q = vertex correction, E v (q = 0) = —E exch . Whereas 
the behavior of the dispersion at short wavelength is dom- 
inated by the exchange self-energy and vertex correction 
terms (see Fig. [3]) , the peak in the dispersion in the long- 
wavelength regime is due to the exchange interaction (the 
RPA term). Furthermore, it is worth pointing out that 
this contribution rapidly increases as one fills more LLs, 
as we will see below. This is a direct consequence of the 
relativistic quantization of the graphene LL spectrum, 
leading to an important LL mixing at higher fillings and, 
as a consequence, building an unusual particle-hole exci- 
tation spectrumP3 

This RPA contribution is absent, however, in the SF 
and SW modes (see Fig. |3j. As a consequence, the 
LL mixing for these modes is less important and makes 
the KH approximation, as the one applied here, a jus- 
tified method (especially at strong magnetic fields and 
for the chemical potential at or near the zero energy 
LL). The results for these modes are also shown in Fig. 
2(a) Electron-electron interactions enter in the disper- 



sion of the former through the term AE^q (q), which 
again due to particle-hole symmetry leads to degen- 
erate modes with contributions Ai^ 1 '"*" 1 / 2 ' ' -1 / 2 ^) = 



E&ijo + Vi d o-A,o(l) and AEi 
V o-x-fi-i{l) respectively. 



-l,-l/2;0, + l/2) 



(?) = £ 



SF 2 



As in the 2DEG, the q -> 
limit of the dispersion of these modes is renormalized 




FIG. 3. (Color online) Decomposition of the ME (full red line) 
and SF (full blue line) mode for v — into the interaction- 
related components, exchange self-energy (dotted line, red for 
ME and blue for SF), vertex correction (dashed-dotted line) 
and RPA term (red dashed line), in units of e 2 /elB- The 
kinetic energy, which yields the same constant offset for both 
modes, is not taken into account in this decomposition. 



makes possible the study of correlation effects by optical 
measurements. 

On the other hand, Larmor's theorem still applies in 
graphene, as one may see from the dispersion of the 
SW mode. This mode has a q — dispersion equal to 
the Zeeman splitting gfigBS z , and a contribution due 
to electron-electron interaction Ai?(°' +1 / 2;0 ^ 1 / 2 ^ (q) = 
^iklo + K?o-o oil)' which is finite only at non-zero wave- 
vectors. This implies that, as in the 2DEG, the g-factor 
is not influenced by the Coulomb interaction. The SW 
disappears if we fill the next LL, for v = 2, with Np = 0, 
v? = 1, and u2 = 1. The dispersions of the ME and SF 
modes (not shown here) are similar to the previous case 
with the difference that the degeneracy of the latter is 
lifted, but only by a constant term equal to the double 
of the Zeeman energy, in agreement with the arguments 
of Sec. HU 



E. Modes at filling v = 4 



The relativistic nature of the LLs in graphene is clearly 
visible if we go beyond Np = 0, as shown in Fig. |4(b) for 
a filling factor of v — 4, with Np = 1, v\ = lj and vl = 0. 
At this filling, the non-equidistancy of the LLs lifts the 
degeneracy of the two ME modes, as well as the two SF 
modes. In addition, the exchange contribution to the 
ME modes, which leads to the peak in their dispersion, 
increases as we decrease the separation between the LLs 
of the electron and the hole. This yields a strong mix- 
ing among the different branches of MEs, as may be seen 
in Fig. |4(a) In fact, the height of the peak associated 
with the MEi [with n e = 2 and = 1, as represented 



from the noninteracting value, EY 



(n e ,n h ) 



g^BBS z . This 



in Fig. 4(b)| , is larger than that of ME 2 (with n e = 1 
and rih = 0). This is due to the linear dispersion of the 
spectrum, which enhances the quantum effects as we go 
to higher filling factors. Taking into account that we are 
showing here only two of the spin-conserving excitations 
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N F =l;v\=l,v\=0 




, = 2 



n = 1- 



SW 



MEi 



SFi 



ME 2 



n = 



SF, 



(b) 

iV F =l;vi=l,v{=0 




FIG. 4. (a) Same as Fig. [2] but for i/ = 4, with 7V F = 1, 
i/| = 1 and = 0. MEi,2 (dashed blue and red lines, respec- 
tively) and SFi,2 (solid blue and red lines, respectively) are 
represented. The thin horizontal lines represent the difference 
in kinetic energy between the electron and the hole E*jb^ and 

^%in - 5 ) Schematic representation of the excitonic modes 
studied. The degeneracy that occurs at Nf = is completely 
lifted at this filling for both, ME and SF modes. For clarity, 
we show in (c) the SF modes separately. 



possible for this filling (the ones involving the more ad- 
jacent LLs to the chemical potential), one can conclude 
that no single MEs will be accessible experimentally at 
finite wave-vectors, but a superposition of them. There- 
fore, the TDHF method in the strong-field approxima- 
tion is not valid for the spin-conserving modes, and the 
inclusion of a much higher number of modes is necessary 



to obtain a reliable result. In fact, this overlap of differ- 
ent MEs leads to a new set of collective modes, the linear 
magneto-plasmons, which have been studied elsewhere.^ 
Notice that the above arguments are valid only at non- 
zero values of the wave- vector, whereas the LL mixing ef- 
fects are less pronounced at q — 0, which is the relevant 
ME energy in magneto-optical experiments!^^ How- 
ever, as we have mentioned above, also at q = the ME 
energy, which is the inter-LL transition energy measured 
in spectroscopy, is renormalized due to electron-electron 
interactions. 

The mixing between different contributions is less dra- 
matic for the SF modes, as one sees in Fig. 4(c) I where we 



show a plot with only the SF modes at this filling. One 
notices how the two modes are clearly decoupled, making 
the use of the KH approximation more justified, because 
of the absence of the RPA term which is responsible, in 
the ME case, for the LL mixing at non-zero values of the 
wave-vector. Although not too clearly, it is appreciable 
that the number of relative extrema (maxima and min- 
ima) in the dispersion of SFi (blue line) is higher than 
for for SF 2 (red line) . This is directly related to the node 
structure of the (hole-) LL wave-functionPSl These max- 
ima and minima lead to hot spots in the dispersion that 
may be detected by Raman spectroscopy techniques!^! 



III. RENORMALIZATION OF THE CHEMICAL 
POTENTIAL 

To gain further insight into the effect of electronic in- 
teractions in a graphene flake, we calculate in this sec- 
tion the exchange correction to the chemical potential, 
from a density-matrix approach. This is the first step 
toward including electron-electron interaction in the sys- 
tem. The correction is intrinsically related to the anti- 
symmetry of the electronic wave-function, which implies, 
even in the absence of interactions, a certain amount of 
correlation between the positions of two particles with 
the same spin. Furthermore, its sign is always negative, 
due to the fact that it is the interaction of each electron 
with the positive charge of its exchange hole. One of 
the effects of Coulomb interaction is a renormalization 
of the chemical potential /x, which at zero temperature 
is the partial derivative of the total energy with respect 
to the number of particles. It contains a contribution 
from the kinetic energy and also from interactions. The 
latter can be written as a mean-field contribution plus 
correlation: // = K + fi ex + fi c , where K is the kinetic 
energy, and fj, ex and \i c are the exchange and correlation 
corrections to the chemical potential, respectively. As 
usual, the direct (Hartree) mean field contribution does 
not appear as it is compensated by the positively charged 
background (or neutralizing background), see the jellium 
model. Furthermore, the exchange interaction can lead 
to a ferromagnetic instability in a dilute electron gasP^ 
In graphene, ferromagnetism due to the exchange inter- 
action between Dirac fermions has also been studied.^ 



In a magnetic field, fi ex can be obtained from the pair 
correlation function g(r) (see Appendix [C| for details of 
the calculation) as 



// / d Y—\g{r) - 1] 
er 



(10) 



where n = 4(1 + N c + N » / '(2irl%) is the electron den- 
sity for graphene in a magnetic field. Np is the in- 
dex of the last occupied LL, related to the filling fac- 
tor by v — ANp + 2, and iV c is a cutoff chosen such 
that (4/V c + 2)N B = 2N U . C _, where N B = A/2ttI% is the 
degeneracy of each LL, A is the surface of the sample, 
N u c is the number of occupied unit cells in the system, 
the factor 2 is due to spin degeneracy, and 4A^ C + 2 is 
the number of filled sub-levels of the valence band for 
undoped graphene. N c is the index of the last LL in 
the band (a kind of bandwidth) and is roughly given 
by N c « N U . C ./(2N B ) = 2nl 2 B /(3V3a 2 ) a 40000/5[T] 
which is always much greater than 1 in practise. The 
fact that N c 1 is just the statement that, with avail- 
able magnetic fields, the flux per unit cell is always much 
smaller than the flux quantum. In this respect, we are 
always in the weak field limit. An exact solution of Eq. 
( 10 ) is possible in the limit N c , Np ^> 1, as shown in Eq. 
( C8 ) . This correction would eventually involve a renor- 



malization of Np, this is, a shift of the chemical potential 
as compared to the non-interacting case. 

Notice that, contrary to the strong magnetic field as- 
sumption done in the previous section, this is the op- 
posite case, namely the weak magnetic field limit. The 
strong field limit is actually the KH app roximation of 
weak LL mixing.^ As stated in Sec. II C the criterion 
for a weak LL mixing is e 2 / (el B co c ) <C 1. In graphene, be- 
cause the fine structure constant ac — e 2 jevp is of order 
one, it means that kpl B <C 1, which means Np ss (i.e. 
Np <C 1) or in other words B 3> 20T. This is the assump- 
tion made in the previous section, whereas in this section 
we assume the opposite limit (B <C 20T or Np ^ 1). In 
a standard 2DEG with a parabolic band, the weak mag- 
netic field limit implies that the typical Coulomb energy 
exceeds the cyclotron frequency up. This allows us to 
start from the Landau Fermi liquid theory at zero mag- 
netic fieldP^l In the case of graphene, this limit is even 
more relevant due to the relativistic quantization of the 
spectrum into non-equidistant LLs, the relative separa- 
tion of which decreases as the energy increases. There- 
fore, even in a strong magnetic field, the strength of the 
Coulomb interaction can be much higher than the sep- 
aration between the LLs adjacent to the chemical po- 
tential (the effective cyclotron frequency in graphene) if 
the system is sufficiently doped. Further simplification is 
possible if 1 <gc Np <§; N c . In this limit we obtain (see 
Appendix [C|) that the exchange correction to the Fermi 
energy behaves asymptotically as 



16y/2 



sin 3n 



(11) 



This contribution is expected since the energy calculated 
above includes the interaction energy of the vacuum of 
negative energy particles. It is interesting to compare 
this leading behavior of the exchange energy, valid for 
high filling factors, to the exchange self-energy obtained 
in Appendix |B| valid at low fillings [see e.g. Eqs. (B8)- 
(B10) for Np = and 1 respectively]. In the two cases 
we obtain the same ~ N\^ 2 leading behaviorPSl Further- 
more, our results agree with the exchange contribution 
calculated for graphene at zero magnetic field, where a 
T, ex ~ —e 2 k c /e contribution was found, k c ~ 1/a being 
an UV cutoff in momenta.^ Taking into account that 
N c ~ (l B /a) 2 , our results for graphene in a magnetic 
field qualitatively agree with those at B = 0. Notice 
that whereas fi ex diverges as N^ 2 for the single parti- 
cle dispersion, the dispersion of a particle-hole pair di- 
verges only logarithmically because the terms propor- 
tional to the square root of N c for each particle cancel 
each other, leading to a behavior E exc h oc \ogN c . This 
divergence can be reabsorbed into a renormalization of 
the Fermi velocity^ and its effect for cyclotron resonance 
measurements has been studied in detail by Shizuya.^l 
This renormalization of the chemical potential due to 
Coulomb interaction should affect the scanning single- 
electron transistor measurements of compressibility in 
grapheneP^ 



IV. SUMMARY AND CONCLUSIONS 

In conclusion, we have studied the SF, SW and ME 
(or magneto-plasmon) modes in graphene in the inte- 
ger QHE regime, in the Kallin-Halperin approximation. 
The ME dispersion in a 2DEG is not renormalized in 
the long-wavelength limit due to Kohn's theorem for sys- 
tems with a parabolic band and Galiean invariance. As a 
consequence, the correction due to the direct interaction 
between the electron and the hole is neutralized by their 
difference in exchange self-energy, E v (q = 0) = —E exc h, 
leading to a dispersion that tends to mwc at zero wave- 
vectorPS In graphene, Kohn's theorem does not apply 
and the dispersion of the ME modes is renormalized due 
to many-body effects even at q = 0. 

On the other hand, virtual transitions from the vac- 
uum (valence band) enhance the depolarization term of 
the spin-conserving ME dispersion, which enters through 
the RPA contribution and which leads to an important 
LL mixing. We have shown that the mixing is higher as 
we increase the LL filling and/or decrease the magnetic 
field, invalidating the applicability of KH approximation 
for v > 2, w hich n eeds to be restricted to the large-field 
Np = case One of our main conclusions is that, for 
ME modes, methods involving more inter-LL transitions 
than only one need to be considered in the calculation of 
the spin-conserving collective excitations. This superpo- 
sition of several inter-LL transitions is at the origin of the 
strongly-dispersing linear magneto-plasmons, which have 
been obtained within an RPA treatment of the electron- 
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electron interactions! 30 * 31 * 

In contrast to the spin-conserving ME modes, the de- 
polarization term is absent in collective excitations where 
the particle and hole components have opposite spin, and 
the amount of mixing is less important. Therefore the 
KH approximation can still be used for these modes in 
undoped or slightly doped graphene in a strong magnetic 
field. In a 2DEG, the zero-wave-vector limit of the KH 
correction of SF modes has a finite contribution, because 
E v (q — 0) = — (l/2)E exc / l in this case. In graphene, the 
dispersion of these modes is also renormalized at zero 
wave-vector and leads to a correction that could be de- 
tected in inelastic light scattering ex perimen ts, by using 
the same techniques as for a 2DEGp ^ 29 * 40 * In contrast 
to Kohn's theorem, we have shown that Larmor's theo- 
rem applies to graphene, so that the q — ¥ limit of the 
SW dispersion is equal to the Zeeman splitting and the 
g-factor is independent of many-body interaction, as in 
a standard 2DEGP3ln addition, the g- factor is also only 
weakly affected by band effects in graphene: the effective 
g-factor was measured to be close to its bare value of 2, 
see Ref. M 

Finally, we have calculated the exchange shift of the 
chemical potential in the weak-magnetic-field limit. We 
have found that, as for strong magnetic fields, the ex- 
change correction to the chemical potential diverges with 

1/2 

the ultraviolet cutoff as ~ N c . However, when the dis- 
persion of an electron-hole pair is considered, the correc- 
tion associated with the difference in exchange self-energy 
between the particle and the hole, diverges only logarith- 
mically. This correction leads to a renormalization of the 
Fermi velocity that seems to explain some recent experi- 
mental resultsPSEH 



Appendix A: Poles of the response function in the 
time-dependent Hartree-Fock and strong field 
approximations 

Within the TDHF approximation, the dispersion rela- 
tion of the excitonic modes is defined by the poles of the 
response function, which are solutions to the eigenvector 
equation [see e. g. Ref. |46]ior the 2DEG] 



{KMDiu)]-* - S s ^S shs *V*s^(ci) 



7,5 



+ S sl!S$ S sf , s ,V^(q)} B 7a (q) = 0,(A1) 



where B 1 g(q) are the basis states, and a = (X a ,n a , s^) 
labels a particle with band index A Q , LL n a , and is 
the z-component of its spin. The sum in Eq. ( |Al| is 
restricted, in the strong-field approximation (i.e. weak 
LL mixing e 2 /(eIb^c) ^ 1), to pairs of indices such that 



n s 



rip = m, and s* 



4 = S * 



This is what we call the Kallin-Halperin approximation. 
The matrix elements of the two-particle propagator are 



D a ,p{u>) 



9lJ>BB(s z p - s z a 

Ml -fa) 



E exch + ig 



g^ B B{s% - s z a ) 



rpexch 
^Pa 



iS 



where f a = (/x — X a VFl] 3 1 y / '2n a ) , <d(x) being the step 
function, and S — >• + . The difference in exchange self- 
energy between the electron and the hole reads 
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E exch = E _ £Q 

where the direct term is^ 



^ A v:'.,. .(0) \ 



7 



V, 



a,P;a'P'\ C l) 



•7(^2)° 



l 

,P-"'.e> ^ ^( A «)^( A /3)^(V)^(A^)«c Al (a), C „(/3);c tI (a'),c,( / 3')(q) 
H,v= 



(A2) 



where d at p- t0t > t p* = 5 nct ,o+^,o+^n a /,o+^n fl ,,Oj b (X) = 1, 
&i(A) = A, c Q (a) — \n a \ and Ci(a) = \n a \ — 1, and 



u a ,p; a 'p'(q) = J drv(r - l 2 B ii z x q)J"*, i/3 , (r)J" Qi/3 (r), 

(A3) 



v(r) = e 2 jer being the Coulomb potential and 



1 



(A4) 



where n < — min(rt a , rip), m = n a — np and e 1 ^ = (x + 
iy)/\x + iy\, and sgn(m) m = 1 for m = 0. The exchange 
matrix elements read 
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V, 



. l/W (q) = -*(V2)<^ 



i 

E 



c^(a),c„(/3);c M (a'),c, 



.(/9')(l). 



(A5) 



where 

V a ,p-a'j3' (q) 



for the MEi mode, whereas E^j 2 = E"/^ 1 as given in 



1 27re 2 



;2 



eq 



^K'^'ih&z x q)-^a,/3(^u 2 x q) 

(A6) 



Appendix B: Exchange self-energy contributions 

In this appendix we give analytical expressions for the 
contributions to AE associated to the difference in ex- 
change self-energy between the electron and the hole. For 



Np 
we obtain 

- el B 
for MEi and 



l,v* = 0, using the notation E 



N c 

E 

n=l 



(4y/n - 3) T (n - 
WV2T(n + 1) 



el i 



(4yn + 3)r(n- \) 
16v^r(n + l) 



(B2) 



for ME2, where N c is a high-energy cutoff. Notice that 
^cku) = ^o-i^o m ^ ne nm it N c — >■ 00. The contributions 
to the spin flip modes are 



J 0;l,0 — ^Oil.O 
„2 



for SFi and 

V S_F 2 

^051,0 _ 



s 



^ME; 

J 0;l,0 
2 



^1;0,0(0) 

(4^ - 3) T (n 



(B3) 



No 

E 

n=l 



16V2r(n + 1) 



VI 



e 
el B 



l,-l;0,o(0) 



(B4) 



7T (4^ + 3)r(n- 1 



for SF 2 . Again, Eg.ft = E 



SFi _ V SF 2 



16-\/2r(ra + 1) 
as the cutoff N r tends to 



infinity. On the other hand, the contribution for the SW 
mode is EqTo = e 2 /(eZs)-\/7r/2, which is cutoff indepen- 



dent. The contributions for Np — 0; 
be expressed in terms of the previously given E 



= 1 can 



MEx- 
0;1,0' 



as E 



ME X 
0:1,1 



^ME 2 
J 0:1,1 



J 0;l,0 



for the ME modes and 



e£*\ = Ef.fj = E*^ 1 for the SF modes. 



Finally, the contributions for Np = = = 0, 



^la.o 



Fig.H 


are 


_ e 2 


1 


el B 


128 




8^ 


71=1 





37V2 - 8J Vtt 7 
\^/n (An - 2^ - 3) + 3] T (n 



128-\/2r(n + 1) 



(B5) 



Eq. (Bl |. For the spin-flip modes we have 



■pSFx 

i;i,o 



B 



3 
4 

E- 

n=l 



(B6) 



fn (An - 2y^n - 3) + 3 T (r 



128^/2 



r(n+l) 



for SFi while the contribution associated to th e sec ond 
mode is Ef.^ 



Sq.^q and coincides with Eq. 

nally, the A^-independent contribution to the Sw mode 

v Np ' ; c V-SW _ e 2 11 nr 
,"f IS iJ 1; i jQ - JJ^- Ye V 2- 

In the following, we calculate the exchange en- 
ergy of the system at low fillings. First, one 
(Bl) notices that the exchange self-energy for undoped 



B3 



Fi- 



1, v+ 



0) 



1,U. 



and simi 
-1 



graphene (Np = 0; v? = 

larly for the filling (Np = 

can be calculated as T, ex — J2n=-N 

-e 2 /(2V2el B ) J2n=i r (" + \) / r ( n + !): which can bc 
summed up exactly to give 



= 1)1 
,(0) = 



s ei = 



el B 2V2 



, £(X+j) 
J r(JV c + 1) 



(B7) 



For N c ^> 1 we obtain the asymptotic behavior 



el B 2^2 



^ + 2WN C + O 



1/2' 



(B8) 



It is useful to express this result by substituting \fJST c by 
its magnetic-field dependence \fW~ c ~ l B /a. By doing so, 
we obtain 



e 

sa 



(B9) 



where # stands for some numerical prefactor and a/l B — 

0. 006 We clearly see that the dominant term, as 
in Eq. (Ill]), is magnetic-field independent. 

A similar result is obtained for doped graphene up to 
the first LL of the conduction band. If (Np ~ l;z/| = 

1, v\ = 0) or (N F = 0;z^ = l,t^ = 1), then E ea; is 
computed as 



£- = ^ ^ 1;n) „(0) 

n=— JV C 



- ^ d 1;n , n (0) + 



E 



l+4^-8nr(ri- |) 
1671 T(n+1) 
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Taking ]V C > 1 we obtain the limiting result 



eh 



1 

7T 



'N r . 



16v/2 



(1.0646 + 7 + In N c ) 



(BIO) 



where 7 is the Euler constant and we have approxi mated 



E„=i (VnT(n - 1/2) /n! - n" 1 ) ~ 1.0646. Eq. flBlpp 
could be accordingly expressed in the form of Eq. ( B9p , 
and the result would be the similar as before: the ead- 
ing term in the exchange contribution to the chemical 
potential does not depend on the magnetic field. 



for the K' (— ) valley, where 



4>n,k(x) 



.e- z *' 2 H n {z). (C3) 



In the previous expression z = x — M b /Ib and H n 
are Hermite polynomial, and we have defined 1* = 
\/(l - S n,o)/2 and 2* = + S n . )/2. One obtains 
therefore 



Appendix C: Correlation function 



The one-particle density matrix for the K valley 
(labeled here by +) can be defined as p+(r, r') — 
Tlia n P+,\,n( r i r ') m terms of the density matrix of 
the n-th LL of the A band of the K valley p+,\, n (r, r') = 

Efc^W^Anfc^). where Hnk( r ) are the ^-valley 
LL wave-function. The wave-function for graphene in a 
magnetic field can be constructed from the correspond- 
ing nonrelativistic LL wave-functions of a 2DEG with a 
parabolic band dispersion. In the Landau gauge, where 
the vector potential is A = (0, Bx,0), they can be writ- 
ten as 



1 _i(y- V ')(x+x') |r-r'| 2 

/<..A.„(r.r'i . ... > "7.- e 



~2wl% 



x l* 2 L„_i 



-'|2 



2/ 2 



-'|2 



2/ 2 



The sum EAiiP+,^>™( r ! r ') m the band and LL indices 
is decomposed into an inter- and an intra-band contribu- 
tions 



Vl 



-iky 



( -i\l* n 4> n -l,k{x) 
K<f>nM x ) 



V 



for the K (+) valley, and 



IT 



-iky 



( 



K4>n,k( X ) 



(Cl) 



(C2) 



E 

n=l 



Nr. 



P+,A=-i,n(r, r') + P+,A=+i,n(r, r'). (C4) 



Furthermore, it can be checked that £„=i lfLl_ x {x) = 

(l/2)L] Vn _ 1 (x), and 2* 2 L° (z) = (1/2)^ (x), 

where iVo = N c ,Np. Therefore, neglecting the Zeeman 
splitting, we obtain for the if -valley one-particle density- 
matrix: 



P+(r,r') = 



2nl 2 B 



rl 



2/ 2 



2/ 2 



+ (7V C _> JV F ) 



(C5) 



where (N c — > iV F ) indicates the replacement of N c 
by ATp. Considering the if'- valley contribution, the 
one-particle density-matrix is obtained as p(r, r') = 
2p + (r,r'). From this, one can obtain the pair correla- 
tion function g(r, r'), which is defined as the normalized 
probability of finding an electron at position r given that, 



at the same time, there is another electron at position r'. 
It can be expressed in terms of the density matrix a^22l 
g(r, r') = 1 — \p(r, r')| 2 /[n(r)n(r')] where n(r) = p(r, r) = 
4(1+N c +Np) / (2irl B ) is the electron density and we have 
used the fact that L°(0) = (ra+a)!/(ra!a!). Setting r' = 
we find 



(C6) 



g(r) 



1 



2e 



2/ 2 



27 2 
zi B 



+ (N c -> N F ) 



where N = 2irl B n(r). By using the asymptotic expres- Bessel function of the first kind, valid for n 3> 1, we 
sion e~ x ' 2 Ll l _ 1 (x) ~ yJn/xJ\ {2^/xn), where J\(x) is a 
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obtain for N F ,N C > 1: 



g(r) ~ 1 



4 



+ i>(N c ,r)+ilj{N F ,r) 



(C7) 



where if)(n,r) = 2nJi (rZ B 1 \/2n) /rl B 1 y^n/2 and we 
have approximated iV c — 1 ~ iV c and iVi? — 1 ~ iVp. 



Using Eq. (C7| into Eq. ( |T0| , with n(r) = n being the 
electron density in the isotropic case, we can obtain an 
expression for the exchange energy per particle in the 
large N c , N F limit with the exact solution 



el b 37r N r 



(N F - N c ) K 



+ (N c + N F ) E 



Nf 

N r . 



+ 4V2 (V c 3/2 + N F /2 



+ &n*> 2 e- N *N c [I (N c ) + h (JVc)] + &Ti 3 ' 2 e~ NF N F [7 (N F ) + h {N F )\ 



2V2' 



.3/2 



(C8) 



where K(ri) and E(n) are elliptic integrals of first and 
second kind, respectively, and I n {z) are the modified 
Bessel functions of the first kind. 
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